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,  .  j  .......  v 

-We-' considers. the  motion  of  a  mixture  of  two  fluids,  with  a  diffusion 

r 

effect  obeying  Pick's  law*  for  the  derivation  of  the  model  (1.1)  see  Section  1 

'■  r  -  •  f  ‘  • 

and  references  [2],  [4], [5]  and  [6].  W*T  consider  ,the  full  non-linear  problem 
-  2 

(i.e.,  we  don't  omit  the  term  A  term  in  equation  ( 1 . 1 ^ ) .  Moreover,  we 

i  i.  ••  -  /»)>■«!  ii  tmtC  i  )*f<- 

don't- assume  that  is  small.  He  prove  the  existence  of  a  (unique)  local 

solution,  the  existence  of  a  global  solution  for  small  data,  and  the 
exponential  decay  to  the  equilibrium  solution*  see  Theorem  A,  Section  1. 
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DIFFUSION  ON  VISCOUS  FLUIDS,  EXISTENCE 
AND  ASYMPTOTIC  PROPERTIES  OF  SOLUTIONS 


Main  Notation 

a 

n  ■  n(x) 

Di'  Dij'  Dt 
I  I ,(  •  ) 

H* 


1  *1 
H„  :  Closure  of  C  (0)  in  H  (0) . 

o  o 

I  lm  :  norm  In  L  (fi). 

2  k  1 

L  ,  H  ,  H0  s  Hilbert  spaces  of  vectors  v  -  *V1 ' v2 ' v3  *  such  that 

2  v  1 

t  L  .v^  e  H  ,vi  e  Hq  (i”1,2,3),  respectively.  Corresponding 
notation  is  used  for  other  spaces  of  vector  fields.  Norms  are  defined 
in  the  natural  way,  and  denoted  by  the  symbols  used  for  the  scalar 
fields. 

Let  us  introduce  the  following  functional  spaces  (see,  for  instance,  [7],  [8]  and 
[12]  for  their  properties)  ” 
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5  {0  K  Hk  :  It 


3n 


0  on  T  and  /  o(x)dx  ■  0  }  ,  k  »  2  . 
-  0 
v  I  {  v  e  C  (fl)  :  div  v  "  o  in  U), 

O 


■  “{vet  t  div  v“0  in  SI,  v»n  ”  0 
T  “  {  v  t  i  div  v  •  0  on  fl). 


n, 


H  and  ▼  are  the  closure  of  v 


2  1 

in  L  (ft)  and  H  ,  respectively.  Moreover 
o 


2  1 

t  “  ■  ♦  G,  where  G  I  {  Vp  i  p  e  H  (J1 )} .  Denoting  by  P  the  orthogonal  projection  of 

2  2 
t  onto  B,  we  define  the  operator  A  =  -PA  on  D(A)  I  ■  nv.  One  has 


(Au,v)  “  (<u,v))  =  l  (D.U  .,D.v  ),  Vu 
i.j  3  13 


e  d( a)  ,  v  e  v. 


The  norms  lol^.  Idol  are  equivalent  in  ,  lolj,  IVAol  are  equivalent  in 

3  2 

and  Ivl^,  1  are  equivalent  in  D(A).  We  define  Ivl^  5  <(v,v))i  the  norms 

lvly,  Ivl^  are  equivalent  in  V. 


L  <0,TjX) 

C(0,T|X) 

U 

1 


v(t,x),  VQ(x) 


p(t,x),  P  (x) 
o 


:  Banach  space  of  strongly  measurable  functions  defined  in  ]0,T[ 

with  values  in  (a  Banach  space)  X,  for  which 

T 

ltl2_  i  /  lz(t)l2  dt  <  +  «. 

L(O.TlX)  0 

:  Banach  space  of  X-vector  valued  continuous  functions  on  1 0, T] 
endowed  with  the  usual  norm  *Z,C(0  y,x) 
s  viscosity  (a  positive  constant). 

:  diffusion  coefficient  (a  positive  constant). 

:  mean-volume  velocity.  Initial  mean-volume  velocity. 

:  density  of  the  mixture.  Initial  density. 

Further, 

m  =  inf  p  (x)  ,  MI  sup  p  (x), 
xefl  °  xS!)  ° 


w(t,x),  p(t,x) 


f(t,x) 


P 


1 

|f!| 


/  p  (x)dx. 

n  ° 


We  assume  that  m  >  0. 

:  pressure.  Modified  pressure 

p  «  it  +  X  v  *Vp  -  X2AP  ♦  X(2u+u’)A  log  p. 
:  external  mass-force. 
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He  denote  by  c,  c,  cQ,  c1#  c2,  •••  positive  constant,  depending  at  most  on  ft  and 
on  the  parasMters  U,  X,  m,  M  and  p.  It  is  easy  to  derive  at  any  stage  of  the  proofs, 
the  explicit  dependence  of  the  constants  on  the  parameters. 

For  convenience  we  sometimes  denote  different  constants  by  the  same  symbol  c. 
Otherwise,  we  utilise  the  symbols  c,  c^,  k  e  *. 

1.  Main  results.  In  this  paper  we  consider  the  motion  of  a  viscous  fluid  consisting  of 

two  components,  say,  saturated  salt  water  and  water.  The  equations  of  the  model  are 

obtained,  for  example,  in  (2],  [4],  [51,  and  [6].  Let  us  give  a  brief  sketch.  Let 

p  1 ,  p2  be  the  characteristic  densities  (constants)  of  the  two  components,  v*^(t,x)  and 
( o ) 

v  (t,x)  their  velocities,  and  e(t,x),  dft,x)  the  mass  and  volume  concentration  of  the 
first  fluid.  He  define  the  density  p(t,x)  =  dp^  ♦  ( 1  — d >  p^,  and  the  mean-volume  and 
mean-mass  velocities  v  S  d  v^1^  +  ( 1  — d)  v^\  w  =  e  v*1^  +  (1-e)v^2\  Then  the 
equations  of  motion  are 

j  p  [D^w  ♦  (w*7)w  -  f]  -pAw  -  (u+u')7  div  w  “  -7w, 

/  div  v-  0  , 

I  D^p  div  (pw)  «  0  . 

On  the  other  hand,  Kick's  diffusion  law  (see  [21)  gives  w  ■  v  -  ip  Vp.  By 
elimination  of  w  in  the  proceeding  equations  one  gets,  after  some  calculations. 


p(Dfcv  +  (v*7)v)  -  llAv  -  X  [(v*7)7p  +  (7p*7)vl  + 

X2  1 

(1.1)  +  -  [ (7p*7)7p  -  -  (7p«7p)7p  +  Ap7p)  »  -  7p  +  pf, 

\  P  P 

DtP  ♦  v7p  -  XAp  “  0  , 
div  v  ■  0  . 

We  want  to  solve  system  (1.1)  in  Qj,  S  ]0,Tl  x  !).  Here  p  is  the  modified 
pressure.  We  add  to  system  (1.1)  the  following  initial  boundary-value  conditions. 


< 

1 

o 

on 

]o,t(  x  r. 

on 

]0,T(  X  r. 

vlt-o  “ 

V  (x) 
o 

in 

a. 

p|t-o  ' 

P  (X) 
o 

in 

n  . 

The  first  two  condition*  mean  that  there  i*  no  flux  through  the  boundary. 

In  (5).  [6]  Kazhikhov  and  Snagulov  consider  the  simplified  system  obtained  from  (1.1) 

2 

by  oeiitting  the  term  containing  X  j  moreover,  they  assume  that 

(1.3)  X  <  . 

Under  these  conditions  Kazhikhov  and  9sagulov  state  the  existence  of  a  local  solution 
in  time  (global  in  the  bidimensional  case). 

In  our  paper  we  take  into  account  the  full  equation  (1.1),  and  omit  the  condition 

(1.3) .  For  this  more  general  case  we  provei  (i)  the  existence  of  a  (unique)  local 

solution  for  arbitrary  initial  data  and  external  face  field)  (ii)  the  existence  of  a 

unique  global  strong  solution,  for  small  initial  data  and  external  force  field.  Moreover, 

if  f  5  0,  the  solution  (p,v)  decays  exponentially  to  the  equilibrium  solution 

(P,0).  More  precisely  we  prove  the  following  result: 

*  2  2  2 

Theorem  K.  Let  v  e  T,p  -  p  e  «...  f  (  l  (0,T»  L  ).  Then  there  exists 

'  -  O  O  N 

t,  e  1 o,t] 

such  that  problem  (1.1),  (1.2)  is  uniquely  solveablc  in  O  .  Moreover 

v  e  l2(0,t1)  h2)  n  c  (o.t^  ▼>,  Dfcv  e  l2  (0,t ^  «>,  p-p  e  l^o.t^h2)  r>  c  (o.t^  hh>. 

DtP  e  L2(0,T1>  H1)  and  m  <  P(t,x)  <  M. 

Moreover,  there  exist  positive  constants  kj,  kj,  and  kj  depending  at  most  on 

*  ( t ) 

0,  u,  X  and  on  the  mean  density  P  such  that  that  if 

(1.4)  'Vi  ♦  IP0-P»2  <  k,. 


* 1  *  Or,  equivalently,  depending  on  the  total  amount  of  mass  1 0 1  p  «  /  p  (x)dx. 
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and 


(1.5) 

L  (0,+-»L*) 

than  the  solution  la  global  in  time.  It  t  =  0,  the  solution  (p,v)  decays 


Ifl  ,  <  k„ 

m  2  2 

L  (0,+-»L  ) 


Iv(t)i,  +  ip(t)-pi,  <  dv  i, tip  -pi  >  e 

I  «  O  I  O  2 


-h3t 


exponentially  to  the  equilibrium  aolutlon  (p,0),  1 ,a . 

(1.6) 

for  avarv  t  >  0. 

Theorem  A  also  holds  for  coefficients  u,X  regularly  dependent  on  p,  v,  provided 

they  are  strictly  positive  and  bounded.  In  a  neighborhood  of  the  range  of  values  of  the 

Initial  data  p  (x),  v  (x>.  This  generalisation  can  be  done  without  any  difficulty, 
o  o 

Moreover,  with  standard  techniques,  one  can  prove  that  the  solutions  have  more  regularity 

m 

(up  to  C  )  if  the  data  are  sufficiently  regular  and  the  usual  compatibility  conditions 
hold. 

2 

Local  existence  in  the  general  case  (i.e.  with  the  X  term  and  without  (1.3))  was 
proved  In  the  lnvlscld  case  by  Beirao-da-Velga,  Seraplonl  and  Valli  In  [1].  a  similar 
result,  in  the  viscous  case  and  for  Si  «  *3,  was  proved  by  Secchi  (11).  for  another 
approach  (concerning  Greffl's  model)  aee  [10]. 

2.  the  linearised  equations.  Wa  start  by  proving  the  following  theorem; Theorem  2. 1  Let 
P(t,x>  be  a  measurable  function  verifying 


(2.1) 


0  <  ■  <  p(t,x)  <  M 


a.e..  In  Qj,,  let  r  e  L  (0,T;H)  and  v  e  V. 


Then  there  exists  a  (unique)  strong  solution  v  of  problem 

r 


(2.2) 


p  o^v  -  pdv  •  -Vp  ♦  r 
div  v  “  0 


v  »  0 
v|t»o 


on  ]0,T[x  T, 


v  (x) 
o 


in  V 
in 

in  0. 


2  2 
Moreover  v  e  L  (0,T|D(A) )  nc(0,T>T) ,  D^v  e  L  (0,T/H)  and 

,2‘3>  M,V,C(0  T;V)  +  "  ,V',22 

C(0,T|V>  LZ(0,T,H) 


♦  lAvl2  «  ulv  I2  +  (~  +  -®r)  IFI2, 

*1  L  (0,T»H)  2M  L  (0,T(H) 
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Proof.  Lot  ua  writ*  aquation  (2.2)  in  the  equivalent  form (2. 4) 

P(pDtv)  ♦  pAv  -  F  ,  V|t-0  “  v0<x>  • 

For  brevity  let  ua  put 

X  -  {vs  v  eL2(0,T»D(») >,  V1  eL2(0,T>  H)  }  . 

Fra*  wall  known  raaulta  (aee  [9],  Vol  I,  Chapter  I  i  Theorem  3.1  with  y  ■  B,  X  - 

D(A),  j  -  Oj  and  (2.42)  Propoaition  2.1)  it  followa  that  X  C*  C  (0,T|  V). 

Me  atart  by  proving  the  a  priori  bound  (2.3 )i  an  eaaential  device  ia  to  Introduce 

parameter  e  in  order  to  conveniently  balance  the  eatinatea.  In  R  take  the  inner 
o 

-1  2 

product  of  (2.4)  with  D^v  +  e^A  v,  £q  >  0.  Since  (v1,  Av)  «  2  D^lvly  one  geta 


(2.5) 


2  u  d  2  2 

>  Id  »r  +  *  ~  Ivi  ,  ♦  e  piAvl  < 
t  2  at  v  o 

<  iFl  ID.  vl  +  £  IFI  I AvI  +  E  M  ID.  vl  lAvl  . 
to  o  t 


By  uaing  the  inequalitiee 
<  4  \  lAvl2  +  v  1  IFl2  and 


Irl  ID^vl  <  4  'm  ID^vl2  +  i 
ID^vIIAvI  <  (4M)  %  I AvI 2  + 


r1  IFl2,  IFl  I  AvI 
-1  2 

U  H  iD^vl  one  gets 


« •  'V2 


+  e  ^  lAvl 
O  2 


.  €  EM 

<  <-  +  — )  IFl2  »  -2 —  ID  vl2 
BP  U  t 


Now  fix  £  •  (4m2)  'up  and  integrate  equation  (2.6)  on  (0,T>.  This  gives  the  a 

o 


priori  bound  (2.3). 

Define  lvlj(  =  "left  hand  side  of  equation  (2.3)",  V  £  L2(0,T;R)  x  ▼,  and 

I (F,v  )l„  I  "right  hand  side  of  (2.3)". 

O  7 

We  solve  (2.4)  by  the  continuity  method.  Define  p®  £  (1-a>  p  +  ap,  a  e  10,1). 


Clearly  PQ  verifies  condition  (2.1),  for  any  a.  Define  Tq  £  (1-a)T  +  aT,  where 

Tv  £  (P (pD^v)  -pAv,  v|t«0>  e  V, 

Tv  £  (P(pDtv)  -pAv,  V|  0  >  e  V. 
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Finally  denote  by  y  the  sat  of  values  a  e  [0,1)  for  which  problem  (2.2)  is  solvable 

in  X  for  avary  pair  (F,v  )  e  V-  Clearly  o  e  Y»  because  for  this  value  of  the 

o 

parameter  aquation  (2.6)  because  the  linearized  Navier-Stokes  equation.  Let  us  verify 
that  Y  is  open  and  closed. 


Y  is  open.  Let  a  e  Y  and  denote  by  G(F,v  )  =  v  the  solution  v  of  problem 
*  o  o 


T  v  “  (F,v  ).  From  (2.3)  one  gets  G  e  LIVtX )  (  )  with  IGI  <  1.  Equation 
a  o  y  fa 

o 

T  t  v  -  (F,v  )  can  be  written  in  the  form 
a  +e  o 

o  _ 

(2.7)  [1-e  G (T-T ) J v  -  G(F,v  ). 

o 

Since  IG (T-T ) I  <  IT-TI  ,  equation  (2.7)  is  solvable  for  |e|  <  lT-TlyV  (by  a 
X#  y  Ay  V  A  #y 

Neuman  expansion) . 

Y  Is  closed.  Let  a  e  Y,  a  ♦  a  ,  and  let  V  be  the  solution  of 
■■  m  m  o  n 

Tv-  (F.v  ).  From  (2.3)  one  has  lv  lu  <  l(F,v  )l ...  Since  X  is  an  Rilbert  space 
a  n  o  nXoy 

n  _ 

there  exists  a  subsequence  ♦  v  e  X»  weakly  in  X.  Prom  T,  T  6  L(Xyy)  one  has 

Tv  ♦  Tv,  Tv  ♦Tv  weakly  in  V.  Hence  T  v  ♦  T  v,  i.e.  T  v  “  (F,v  ). 
v  v  a  v  a  a  o 

vo  o 

Let  us  now  return  to  problem  (1.1).  Define 

(2.8)  F(p,v)  S  P  {-p(vV)v  +  X  [(vV)Vp  +  (7p*V)v)  + 

A2  1 

+  —  [ (Vp« V)Vp-  f  ( Vp*7p)7p  +  Ap7p)  +  pf)  . 

P  P 


For  convenience  we  will  use  in  the  sequel  the  translation 

(2.9)  P  -  P  +  o. 


Recall  that  P  Is  a  given  constant.  To  solve  problem  (1.1)  and  (1.2)  in  our 

2  2 

functional  framework  is  equivalent  to  finding  v  e  L  (0,T|D(A)),  vi  e  L  (0,TjH)  and 

2  3  2  1 

a  e  L  (0,T|  H^),  O'  e  L  (O.Ti  H  )  such  that 

P(pD^v)“pAv  -  F(p,v), 


(2.10) 


I  t-0 


v  (x) , 
o 


D^o  -XVo  “  -v*Vo, 


1 1-0 


o  (x), 
o 


(2)  Banach  space  of  linear  continuous  operator  from  V  into  X,  with  norm  I  I 


V.X- 


-7- 


are  given.  Note  that  from  the  above  conditions  on 


where  vq  e  V  and  0q(x)  “P  e  H^lfl), 
c  it  follows  that  a  e  C(0,T»  H2). 

We  solve  (2.10)  by  considering  the  linearized  problem 


P(pD^v)-uAv  -  P< p,v)  5  F, 


(2.11) 


|t-0 


»  v  (x). 


D^O  -  XA<J  “  -v*Vo.  i 

°|t-0  “  °o<x)' 

and  by  proving  the  existence  of  a  fixed  point  (p,v)  «  (p,v)  for  the  map  (p,v)  ♦  (p,v) 
defined  by  (2.11). 

In  order  to  get  a  sufficiently  strong  estimate  for  the  linearized  equation  (2. 11  )j  we 
take  in  account  the  particular  form  of  the  date  v Vo.  As  for  estimate  (2.3)  we  will 
introduce  a  balance  parameter  c  >  0. 

Theorem  2,2  Assume  that  v  e  L2(0,Tii2)nc(0,TiIj)  and  that 
a  e  L2(0,T>h^) nc(O.TjH^).  Then  the  solution  o  e  L2(0,T»H3),  o'  e  L2(0,T»H1)  of 
problem  (2.11)^,  (2.11)^  verifies  the  estimate 

<2*12>  ,O'?(0  T,  R2)  *  ,0,22  3  ‘ 

C<0.T,  Hjj  L2(0,T|H2> 


(  c Jo  1^  +  C  e"3T[lvl6 


+  IVol6 


C(  0, t;H  > 


,  1 

C( 0, T;H  ) 


+  C3  e  [ivl2 

L  (0,T:H  > 

for  every  positive  e  verifying 
(2.13) 


+  >’P*22  2  ]. 

L  (0,T>H  ) 


.  * 

c  <  — 
2c 


where  c^  is  the  constant  in  (2.16).  Here  c^,  Cj,  Cj  are  positive  constants  depending 
only  on  fl. 

Proof.  The  existence  of  a  solution  a  in  the  required  space  follows  with  standard 
techniques  from  the  a  priori  bound  (2.12),  or  from  [9],  Volume  II,  Chapter  4,  Theorem  5.2, 
with  H  ■  H  .  Let  us  prove  (2.12): 
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By  application  of  the  operator  A  to  both  aidea  of  (2. 1 1  >3,  by  multiplication  by 
Acr#  and  by  integration  over  SI  one  gets  V2  Dt  lAol2  +  ( V( XAo-v*Vo ) ,  VAc)  -  0.  Hence 

(2.14)  4  -TT  >Ao|2  +  XI7AOI2  < 

2  dt 

<  (IDvDol  +  IvD2cI)  I7AoI  . 

16  " 

By  using  Sobolev's  embeding  theorem  H  C*L  and  Holder's  inequality  it  follows  that 

( 2. 15)  7  —  IAol2+  XIVAol 2  < 

2  at 

—  V*  —  1A  — 

<  CdDvl  2  iDvl'2  I7al  1  + 

-1/  _  1/ 

+  lvl1  IVol'2  IVAal  2  )  1 VAol  . 

-3  4  2  4 

An  utilisation  of  abc  <  (8e  )a  ♦  (e/2)  b  +  (e/2)  c  ,  e  >  0,  leads  to 

(2.16)  4  lAol 2  +  XI VAol 2  <  c  ElDvl2  + 

2  dt  i 

♦  c1  ElVAol2  +  -j  I*}  IVol2  + 

+  c  e  I  VAol 2  Sj  Ivl2  IVol  *  . 

Hence  for  e  verifying  (2.13)  one  has 

(2.17)  lAol2  +  X  l7Aol 2  <  -§  (Ivl2  I  Volt  + 

dt  e 3  i  i 

♦  Ivl*  IVol2)  +  c9C  ( I  Dvl  2  +  I  VAol 2  ) , 

where  the  constants  c  depend  only  on  SI.  This  proves  inequality  (2.12).  Recall  that 
lol2  <  clAol,  lol3  <  cIVAol  . 

Remark.  One  could  also  consider  the  linearized  equation  o  +  v*Vo  -  XAo  -  0  instead 
of  (2.11>9>  then  estimate  (2.17)  holds  with  o  replaced  by  o  and  without  the  term 
c^e  IVAol  .  In  this  case  the  solution  o  of  the  linearized  problem  verifies  the  maximum 
principle  (which  doesn't  hold  for  the  solution  of  <2.11)3).  However,  the  linearization 
(2. 1 1 ) j  seems  to  be  more  in  keeping  with  the  linearization  (2.11)^  Besides,  the  maximum 
principle  will  be  recovered  for  the  solution  of  the  full  nonlinear  problem  (2.10)3. 

3.  The  nonlinear  problem.  Local  existence.  We  will  not  take  care  of  the  explicit 
dependence  on  u,X,m,M>  some  of  the  constants  c,ck,  depend  on  these  fixed  quantities.  In 


order  to  simplify  the  equations,  we  denote  by  Kq,  K^, 


constants  depending  on  the 


In  thin  section  we  solve  (2.10)  by  proving  the  existence  of  a  fixed  point 


(P,v)  -  (p,v)  for  system  (2.11).  Define 
K  i  {v  :  v 


1 


|t-0  _  V0(X)'  '"'V  T  2 

L  ( 0 ,  T  j  H  ) 


♦  Ivl2  «■  Iv’l2  <  2  C.  Iv  I2  } 

C(0'T,T)  L2(0,T,H)  4  oV 


*2  1  {0  !  °lt-0  _  °o(x)' 


L‘(0.T»H')  C(0,T|Hk) 


<  2  c,  lo  if.  lo'l  ,  <  K  , 

2  0  2  ^(O.T.H1)  ° 


lo  -  o  I  <  r  )  , 

“c^)  2 


where  C.  =  lifmin  { p ,  M ,  (  bxj 2 ) /  ( 4M2 ) }  ]  1  and  K  =  \/2c  lo  I  +  C  / 4c  c  Iv  I  lo  I  . 

4  L  1  O  2  o  2  24oVo2 

Here  c  -  c(fl)  is  a  positive  constant  such  that 

(3.1)  Iv-wl,  <  c  lvlv  lwl2  ,  V  V  e  v,  w  e  h2  . 

Note  that  for  every  n  6  Ij  one  has  in 

(3.2)  J  <  P(t,x)  <  M  +  ~  . 

we  now  evaluate  the  I.2  norm  of  F  S  F(p,v).  By  using  Sobolev's  embeding  theorem 
H'ceL6  end  Holder's  inequality  one  easily  gets 


(3.3) 


lF(P,v)l2  <  C  Ivl2  I  Dvl  1  +  C  Ivl2  ID2ol  • 

•  lD2CJl  ,  +  C  Ivl,  I  Dvl  ,  I  DCl2  + 
+  C  I  Do  I  2  I  D^ol  1  +  C  IDOl®  + 
+  C  I  f  1 2  . 


Consequently, 
(3.4)  lF(p,v)l2. 


L  (0,T|H) 

♦  C  lo  1^  T  +  C  •  f  1 2 

°  2  L2(0,T)H) 


<  C  dv  I  +lo  l-)7/2  TV2 
o  V  o  2 


Hence*  by  using  (2# 3),  it  follows  that  the  solution  v  of  (2»11)j*  (2*  11  >2  verifies 
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(3'5>  ,V,C(0,T.V)  +  ,  2  ♦  '"‘‘V  ,  2 

L  (0,TjH  )  L  (0#T;Xt  > 


<  C.  Iv  1 3  +  K  (/t  +  T)  ♦  C  Ifl2 

4  °  v  1  l2<o.t,«>  . 


On  the  other  hand,  from  (2.12), 

)  lol2  +  lol2 

C<0,T|Hn>  L  (0,T(Hjj> 


<  {  C,  Id  I?  +  K,  e“3T  +  K.  E  } 
t  O  i  3  4 

-1  2 

Now  we  fix  E  >  0  such  that  K  E  <  2  C,  Id  I,  .  I 

4  2  o  2 


Finally,  by  fixing  a 


sufficiently  small  T  >  0,  it  follows  that  v  e  ^ ,  o  t  Ij.  The  estimate  for  0  o 
follows  by  using  equation  (2.11 >3  and  (3.1).  The  estimate  for  the  sup  norm  of 
d  -  in  Qt  is  proved  as  follows: 

Clearly,  lo(t)  -  o  I  <  J  lo'(s)l,  ds  <  K  T1/2  .  On  the  other  hand  lo(t)  -  o  I 

Olio  o  T 

1/3  °  2/3  C(0) 

<  Cj  loltj-o^^  lo(t)-Ool2  ,  where  Cj  depends  only  on  fif  recall  that 

(fl)  C^C(Si).  Consequently 

lo  -  a  I  <  C,  r1/3  T1/6  <c,  /3cT  lo  I,  *  I a  l,)2/3  , 

°  C(Qt>  5°  6  2  0  2  0  2 

2 

where  C  "  C,(U)  is  a  positive  constant,  such  that  lol  <  C,  lol,  ,  Vo  e  H„(fl). 

6  6  C(0>  6  2  M 

Hence,  by  choosing  (if  necessary)  a  smaller  value  for  T,  one  gets  lo  -  o  I  <  m/2  . 

°  CtQ^,) 

Now  we  utilize  Shauder'a  fixed  point  theorem.  Clearly  K  i  ♦  Kj  is  e  convex, 

2  2  2  2  —  — 
compact  set  in  I.  (0,T;I.  )  x  L  (0,T;L  ).  Let  us  denote  by  ♦  the  map  4(p,v)  -  (p,v), 

defined  by  (2.11).  Since  *  (K)  C  K,  it  is  sufficient  to  prove  that  *  t  K  ♦  K  is 

2  —  —  2  —  ~  2 

continuous  in  the  L  topology.  If  vn  ♦  v  in  L  (Q^),  Pn  ♦  p  in  L  (QT>,  is  follows  by 

-  -  2  2  12 
compactness  arguments  that  v  v  weakly  in  L  (0,T;H  )  and  in  H  (0,T|I.  ),  and  that 

n 

—  —  2  2  12  — 

Vp  ♦  Vp  weakly  in  L  (0,T:H  )  and  in  H  (0,TiL  ),  In  particular  P  is  a  bounded 

V2«,  2-E2  a-  3  " 

sequence  in  H  (0,T;H  >  C^C0'01^}^)  *  *,  for  suitable  positive  E^,  a. 


(  )  o-Holder  continuous  functions  in  Q 
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Hence  ♦  p  uniformly  in  Moreover,  v^  and  Vp^  are  bounded  in  H^tO.T/l1)  and 

in  H/^(0,T)H%  respectively.  Hence  v  ♦  v  and  Vp  ♦  Vp  strongly  in  the  L* 

n  n 

topology.  It  follows  from  (2.8)  at  F(p  ,  v  )  ♦  F(p,v)  as  a  distribution  in  Q_j 

n  n  * 

consequently  F(p  ,  v  )  F(p,v)  weakly  in  L2(Q  ),  because  F(p  ,v  >  is  a  bounded 
n  n  t  n  n 

sequence  in  this  space.  Analogously,  v  •  Vp^  ♦  v  •  Vp  strongly  in  L^Q^).  It  follows 

2  2 

from  the  linear  equations  (2.11)  that  v  ♦  v  and  p  *  p  in  L  (0  )  and  L  (Q_ ) , 

n  n  T  T 

respectively.  Hence  *  is  continuous.  This  finishes  ths  proof  of  the  existence  of  a 
local  solution.  Uniqueness  will  be  proved  in  Section  5. 


4.  Global  solutions.  Asymptotic  behavior. 

In  this  section  the  constants  Cg  depend  at  most  on  (2  and  on  the  quantities 
U,  X  and  p,  i.e.  on  the  total  amount  of  mass  |Qj  p.  He  assume  that 

(4.1)  lo  I,  <  (2C)_1  p, 

O  4 

hence  p/2  <  m  <  M  <  3  p/2.  Let  (P,v)  be  a  solution  of  (1.1).  Proa  (2.6)  for 
2  -1 

e  *  (4M  )  my,  and  from  (3.3)  one  gets 
o 


(4.2) 


2  '  V 


2  dt 


V  '  2 

8  M 


lAvl2  < 


<  C  (Ivl2  +  loi2  )  (lvl2  +  lolj)  +  c  lol 2  ♦  c  Ifl2 


-1, 


where  C  depends  only  on  0,  u,  p.  On  the  other  hand,  from  (2.17)  for  e  ■  ( 2Cg )  X, 


one  gets 

(4.3)  lAol2  +  |  IVAol2  <  C  (Ivl®  +  lol®  )  . 

From  (4.2)  and  (4.3)  it  easily  follows  that 
f  *vl2  ♦  lAol2)  +  §  lDtvl2  + 

2 

+  S11-—  I AvI 2  +  j  IVAol2  <  C(lvl®  +IAol6+lf|2)  . 
16M 
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(4.4) 


In  fact,  if  Cn  (  lv(t)ly  +  IAd<t>2!)2  -  C10/2  it  must  be,  from  (4.4),  that 
“  dv(t)l2  +  I Ad(t ) I2 )  <  0. 

at  V 

Let  us  now  prove  the  last  aasertion  in  theorem  A.  Under  the  hypothesie  (4. 5) ^ ,  it 
follows  from  (4.4)  that 

■g  (Ivl2  ♦  I  Adi 2 )  <  -  (Ivl*  +  I  Adi2)  . 

this  proves  (1.6)  .  q 

5.  uniqueness .  We  prove  that  the  solution  (0,v)  of  problem  (1.1)  is  unique  in  the  class 
in  which  existence  was  proved;  see  Theorem  A.  We  remark  that  more  careful  calculations 
lead  to  uniqueness  in  a  larger  class. 

Let  (P,v>,  (p,v)  be  two  solutions  of  problem  (1.1),  (1.2)  and  put 
u  ■  v-v,  n  ■  p-p.  By  subtracting  the  equations  (2.10)^  written  for  (p,v)  and  (p,v> 
respectively,  and  by  taking  the  inner  product  with  u  in  H  one  gets 

j  ■—  (Pu.u)  +  U  lul2  “  _  j  (v*Vp,u2>  + 

+  j  (Ap,u2)  -  (y.D^v'u)  +  (F-F,u>  . 

2  2 

By  using  (Ap,u  )  “  -(Vp,Du  >  ,  we  show  that 
(5.1)  ~  (pu.u)  +  |  lul2  <  J  Ivl^lVpIJul2  * 

♦  ^  I7pl 2  lul2  +  ClD^vl 2  lul2  + 

+  j  I Anl2  +  (r-F,u)  . 
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On  the  other  hand,  by  subtracting  aquation*  ( 2 . 1 0 )  3  written  (or 

—  —  2 
(P,v)  and  (or  (p,v)  raapactivaly,  and  by  taking  tha  innar  product  with  An  in  L  (Q) 

ona  gata 

(5.2)  4  '’H2  *  T  «&nl 2  <  C  IVpl2  lul2  ♦ 

+  c  ivi2  ivni2  . 

By  adding  (5.1)  and  (5.2)  it  (ollowa  that 

(5.3)  “  t(pu,u)  *  IVnl2]  +  li  lul2  +  j  lAnl2  < 

<  e  (t)  (lul2  ♦  ivni2)  +  (f-f,u), 
whara  9,(t)  ia  a  raal  integrabla  (unction  on  10, TJ. 

« 

On  tha  othar  hand,  by  uaing  Sobolev ‘a  eabeeding  thaoreas  and  Holder' a  inequality 
2-12 

(and  alao  ab  <  ea  +e  b  )  the  reader  aaaily  veri(iea  that  given  t  >  0  there  exist*  an 
integrable  real  (unction  6^ (t)  (dependent  on  p,  p,  v,  v  and  on  e)  such  that 

(5.4)  |  (F— F,u)  |  <  62<t)  lul2  ♦  e  (lnl2tlul2  )  . 

2  -1 

By  using  lul  <  m  (pu,u)  ,  (5.3)  and  (5.4)  it  follows  that 

^  ( ( pu , u)  +  IVnl2]  <  (0,(0  +■  e2(t)>  ( ( pu , u )  ♦  IVnl2)]. 

Uniqueness  Follows  now  (roe  Gronwall's  lemma  and  (rom  U|t_Q  “  0,  n | “0. 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  we  consider  the  Motion  of  a  continuous  medium  consisting  of 
two  components,  say  water  and  a  dissolved  salt,  with  a  diffusion  effect 
obeying  Pick's  law.  We  denote  by  v,w,p ,w ,u,X  the  mean-volume  velocity,  mean- 
mass  velocity,  density,  pressure,  viscosity  and  diffusion  constant, 
respectively.  By  using  Pick's  law  we  eliminate  w  from  the  equations  and  we 
obtain  (1.1),  where  p  is  the  modified  pressure;  see  Section  1  and  references 
[2],  [4], [5]  and  [6].  The  initial  boundary  conditions  are  given  by  equation 
(1.2). 

Kazhikhov  and  Smagulov  [5], [6]  consider  equation  (1.1)  for  a  small 

diffusion  coefficient  X.  More  precisely,  they  assume  that  condition  (1.3) 

2 

holds;  moreover,  they  omit  the  X  term  in  equation  (1.1),,.  Older  these 
conditions  they  prove  the  existence  of  a  unique  local  solution  for  the 
3-dimensional  motion  (in  the  bi-diaensional  case,  solutions  are  global). 

In  our  paper  we  consider  the  full  equation  (1.1),  without  assumption 
(1.3),  and  we  provei  (i)  the  existence  of  a  (unique)  local  solution;  (ii)  the 
existence  of  a  global  solution  in  time  for  small  initial  velocities  and 
external  forces,  and  for  initial  densities  near-constant;  (iii)  the 
exponential  decay  (when  t  ■*■  +•)  of  the  solution  (p,v)  to  the  equilibrium 

A 

solution  (p,0) ,  if  f  =  0.  See  Theorem  A,  Section  1. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


DIFFUSION  ON  VISCOUS  FLUIDS,  EXISTENCE 
AND  ASYMPTOTIC  PROPERTIES  OF  SOLUTIONS 


Main  Notation 

a 

n  ■  n(x) 

V  Dij'  Dt 

H* 


H.  Beirao-da-Veiga* 


t  an  open  bounded  aet  In  Br,  locally  situated  on  one  side  of  Its 
boundary  T,  a  regular  (say  C*)  manifold. 

>  unit  outward  normal  to  T. 

«  3/3x1,  32/3xi3xj,  3/3t. 

t  norm  and  scalar  product  in  L2(ft). 
k  2 

i  Sobolev  space  H  '  (fl)  with  norm 

lol2  S  l  ID1^2, 

1-0 


where 


Further, 


lD18l2s 


l 

|«|«1 


IDaol2, 


I D^Ol 2  s  l  ID°Ol2  . 
*  l«|-l  “ 


*0 

I  l_ 

L2, 


t  Closure  of  C  (0)  in  H  (Q). 
o 


I  norm  In  L  (0). 

j  Hilbert  spaces  of  vectors  v  «  (v1 ,v2 , v3 )  such  that 

2  Jr  1 

Vj  (  l  ,vi  e  H  ,v^  e  Hq  (1-1,2,3),  respectively.  Corresponding 
notation  is  used  for  other  spaces  of  vector  fields.  Norms  are  defined 
in  the  natural  way,  and  denoted  by  the  symbols  used  for  the  scalar 
fields. 

Let  us  Introduce  the  following  functional  spaces  (see,  for  instance,  [7],  [8]  and 
(121  for  their  properties)  « 
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H*  S  {o  K  Hk  «  f2  -  0  on  r  «nd  /  o  ( x ) dx  -  0  >  ,  k  >  2 

"  3n  „  n 

v  s  {  v  e  C  (fl)  i  div  v  “  o  in  0}, 

O 

2 

■  -{vet  :  div  v  “  0  in  Q,  v*n  “  0  on  I*}, 


■  -  {  v  e  Hg  i  div  v 


0  on  0} . 


■  and  V  are  the  eloeure  of  v 


2  1 

in  L  (fl)  and  H  ,  respectively.  Moreover 
o 


2  1 

L  -  H  +  G,  where  G  5  {  Vp  :  p  e  H  (0)}.  Denoting  by  P  the  orthogonal  projection  of 

2  2 
t>  onto  H,  we  define  the  operator  A  s  -PA  on  D(A)  2  ■  nv.  One  has 


Vue  D(A),  v  e  V. 

2 


( Au / v )  -  ((u,v))  2  l  (D  0  ,D  v  ), 

i.j  3  3 

The  norma  lolj,  Idol  are  equivalent  in  H^,  lolj,  IVAol  are  equivalent  in 
and  Ivlj,  lAvl  are  equivalent  in  D(A).  We  define  Ivl^  =  (<v,v))»  the  norms 


Ivly,  Ivl 1  are  equivalent  in  V. 


L  <0,T»X) 

C<0.T»X) 

U 

1 


v(t,x),  vQ(x) 

p(t,x),  P  (x) 
o 


:  Banach  space  of  strongly  measurable  functions  defined  in  ]0,T[ 

with  values  in  (a  Banach  space)  X,  for  which 

T 

Izl2  5  /  '*(t)|2  dt  <  + 

L (O.TlX)  0 

x  Banach  space  of  X-vector  valued  continuous  functions  on  [0,T] 
endowed  with  the  usual  norm  lzlc((J  r,X). 
t  viscosity  (a  positive  constant). 

:  diffusion  coefficient  (a  positive  constant), 
i  mean-volume  velocity.  Initial  mean-volume  velocity, 
s  density  of  the  mixture.  Initial  density. 

Further, 

m  S  inf  p  (x)  ,  M2  sup  P„(x), 

xen  °  xen  ° 


*(t,x),  p(t,x) 


f <t,x) 


P 


1 

mi 


/  p  (x)dx. 

n 


We  assume  that  m  >  0. 
s  pressure.  Modified  pressure 

p  «  v  +  X  v  •  Vp  -  l2Ap  ♦  X(2if+y')A  log  p. 
t  external  mass-force. 
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N*  denote  by  c,  e,  cQ,  c,,  Cj,  •••  positive  constant,  depending  at  most  on  n  and 
on  the  paraawtara  U,  X,  a,  M  and  p.  It  la  aaay  to  darlva  at  any  stags  of  the  proofs, 
ths  explicit  dependence  of  the  constants  on  the  parameters. 

For  convenience  we  soaetlaes  denote  different  constants  by  the  saae  symbol  c. 
Otherwise,  we  utilise  the  symbols  c,  c^,  k  e 

1.  Main  results.  In  this  paper  we  consider  the  motion  of  a  viscous  fluid  consisting  of 
two  components,  say,  saturated  salt  water  and  water.  The  equations  of  the  model  are 
obtained,  for  example.  In  (2],  [4],  (5],  and  [6].  Let  us  give  a  brief  sketch.  Let 
p  1 ,  p2  be  the  characteristic  densities  (constants)  of  the  two  components,  v^(t,x)  and 
v(2,(t,x>  the^r  velocities,  and  e(t,x),  d(t,x)  the  mass  and  volume  concentration  of  the 
first  fluid.  He  define  the  density  p(t,x)  =  dp 1  +  (1-d)  ,  and  the  mean-volume  and 

mean-mass  velocities  v  =  d  v  +  (1-d)  v  ,  w  =  e  v  +  (1-e)v  .  Then  the 

equations  of  motion  are 

C pfD^w  +  (w*7)w  -  f)  -U& w  -  (p+p')7  dlv  w  “  -7», 

/  dlv  v-  0  , 

I  Dfcp  ♦  dlv  (pw)  »  0  . 

On  the  other  hand,  Tick's  diffusion  law  (see  [2))  gives  w  *  v  -  Xp  1Vp.  By 
elimination  of  w  in  the  proceeding  equations  one  gets,  after  some  calculations, 

p(Dfcv  +  (v*7)v)  -  Uiv  -  X  l(v*7)7p  +  (7p*7)v)  + 

X2  1 

(1.1)  +  -  ( (7p*7 )7p  -  —  (7p»7p)7p  +  Ap7p]  -  -  7p  +  pf, 

\  P  P 

Dfcp  +  vVp  -  XAp  “  0  , 
dlv  v  ■  0  . 

He  want  to  solve  system  (1.1)  in  Qj  =  JO.Tt  x  0.  Here  p  Is  the  modified 
pressure.  He  add  to  system  (1.1)  the  following  initial  boundary-value  conditions. 
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v  “  0 


on  ]0,T(  x  f 


|£-  o 

3n 

on 

]o,t[  x  r. 

v|t-0  ’  Vo(x) 

in 

fl. 

P|t-0  "  Po(K> 

in 

n  . 

(1.2) 


The  first  two  conditions  mean  that  there  is  no  flux  through  the  boundary. 

In  (51.  [6]  Kaxhikhov  and  Smagulov  consider  the  simplified  system  obtained  from  (1.1) 
by  omitting  the  term  containing  12>  moreover,  they  assume  that 
d.3>  X<£_. 

Under  these  conditions  Kashikhov  and  Smagulov  state  the  existence  of  a  local  solution 
in  time  (global  in  the  bidimensional  case). 

In  our  paper  we  take  into  account  the  full  equation  (1.1),  and  omit  the  condition 
(1.3).  For  this  more  general  case  we  provei  (i)  the  existence  of  a  (unique)  local 
solution  for  arbitrary  initial  data  and  external  face  field*  (ii)  the  existence  of  a 
unique  global  strong  solution,  for  small  initial  data  and  external  force  field.  Moreover, 
if  f  =  0,  the  solution  (P,v)  decays  exponentially  to  the  equilibrium  solution 

(p,0).  More  precisely  we  prove  the  following  result: 

*  2  2  2 

Theorem  h.  Let  v  e  V,p  -  p  e  K.,  f  61  (0,T>  L  ).  Then  there  exists 

— — — — ’  O  O  w  '  -  -  - 

t1  e  ]0. Tl 


such  that  problem  (1.1),  (1.2)  is  uniquely  solveable  in  •  Moreover 

2  2  2  “^2  3  2 

v  e  l  (o,Tt j  ■  )  n  c  (o,t^ :  v),  Dfcv  e  l  (O.t^  ■>,  p-p  e  l  <0.'Vh,i>  n  c  (o.t^  hn>. 

0tP  e  L2(0,T1i  H1)  and  m  t  p(t,x)  <  M. 

Moreover,  there  exist  positive  constants  k1#  k2,  and  kj  depending  at  most  on 

*  ( i ) 

0,  u,  1  and  on  the  mean  density  p  such  that  that  if 

(1.4)  lv  I.  ♦  Ip  -pi  *  kt, 

o  1  o  2  1 


(1)  Or,  equivalently,  depending  on  the  total  amount  of  mass  I  Cl  t  p  »  /  pQ(x)dx. 
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(1.5)  Ifl  „  ,  <  k  , 

L  (0,+**jL  ) 

than  the  solution  la  global  in  time.  If  f  =  0.  the  aolutlon  (p,v)  decays 

exponentially  to  the  equilibrium  solution  (p,0),  l.e. 

-k,t 

(1.6)  !▼(«!,  +  lp(t)-pl2  <  (,v<J11+IP0-pl2>  «  . 

for  every  t  >  0. 

Theorem  A  also  holds  for  coefficients  U,X  regularly  dependent  on  p,  v,  provided 

they  are  strictly  positive  and  bounded,  in  a  neighborhood  of  the  range  of  values  of  the 

initial  data  p  (x),  v  (x).  This  generalisation  can  be  done  without  any  difficulty, 
o  o 

Moreover,  with  standard  techniques,  one  can  prove  that  the  solutions  have  more  regularity 
(up  to  C**)  if  the  data  are  sufficiently  regular  and  the  usual  compatibility  conditions 


Local  existence  in  the  general  case  (i.e.  with  the  X  term  and  without  (1.3))  was 
proved  in  the  inviscld  case  by  Beirao-da-Veiga,  Serapionl  and  Valli  in  [1].  a  similar 
result,  in  the  viscous  case  and  for  Q  »  was  proved  by  Secchi  [11].  For  another 

approach  (concerning  Greffi's  model)  see  [10]. 

2.  The  linearised  equations.  Me  start  by  proving  the  following  theorem : Theorem  2. 1  Let 
p(t.x)  be  a  measurable  function  verifying 

(2.1)  0  <  a  <  p(t,x>  <  M  ,  a.e. .  in  Op,  let  F  e  L2(0,TiM)  and  vq  e  ▼. 

Then  there  exists  a  (unique)  strong  solution  v  of  problem 


p  D^v  -  v4v  »  -Vp  ♦  F  in 
dlv  v  ■  0  in  C^, 


on  ]0,T(x  T, 


V|.  „  “  v  (x) 
lt-0  o 


2  2 
Moreover  v  e  L  (0,T»D(A) )  nc(0,T»V),  Dtv  e  L  (0,T»M>  and 


M,V,C(0  Ti V)  +  ■  ,V’,22 

C(0,T»V)  L  (0,T»H) 


♦  lAvl2,  <  U,Vo'v  +  +  ~&2)  ,F,22 

4M2  L2(0,T>H)  2M2  LZ(0,TiH) 


Proof.  Lot  us  writ*  aquation  (2.2)  in  th*  equivalent  form(2.4) 

P(PDt*)  +  WAV  -  F  ,  V|fc_0  -  vq(X)  . 

For  brevity  l*t  us  put 

X  -  {v:  V  «L2(0,T|D(A))#  V*  eL2(0,T»  H)  }  . 

Pros  well  known  result*  (•••  [9],  Vol  I,  Chapter  I  i  Theorem  3.1  with  y  "  H,  X  * 

D(A) ,  j  -  Oi  and  (2.42)  Proposition  2.1)  it  follows  that  X  Cj  C  (0,T»  V). 

W*  start  by  proving  the  a  priori  bound  (2.3)t  an  essential  device  is  to  Introduce  a 

parameter  e  in  order  to  conveniently  balance  the  estimates.  In  H  take  the  inner 

-1  2 

product  of  (2.4)  with  Dfcv  +  e^A  v,  t,Q  >  0.  Since  (v',  Av)  «  2  Dt,v,v  one  9®tB 

(2.S)  "  'V'2  +  2  £  lvtl  +  Eo'J,AV'2  * 

<  IfI  ID  vl  ♦  e  Ipl  lAvl  +  £  M  ID  vl  lAvl  . 
to  o  t 


By  using  the  inequalities  IFI  IDtvl  <  4-1m  lDtvl2  +  m  1  IFI2,  IFI  lAvl 
<  4-1w  lAvl2  ♦  W-1  IFI2  and  iD^llAvl  <  (4M)_1W  lAvl2  +  u_1M  IDtvl2  one  gets 


(2.6) 


i  ■  IDtvl2  ♦  ^  — 


2  dt  +  %  2  ‘AV*2  « 


,  e  cM  , 

« «i ♦  f >  ,FI  ♦  -V  'V  • 

2  *1 

How  fix  e  -  (4m  )  n u  and  integrate  equation  (2.6)  on  (0,T).  This  gives  the  a 
o 

priori  bound  (2.3). 

2  2 
Define  lvlj(  =  "left  hand  aide  of  equation  (2.3)",  V  s  L  (0,t>M)  x  ▼,  and 

I (F,v  )lj  =  "right  hand  side  of  (2.3)". 

we  solve  (2.4)  by  the  continuity  method.  Define  pa  r  (1-a)  p  +  ap,  a  e  (0,1). 

Clearly  Pa  verifies  condition  (2.1),  for  any  a.  Define  Ta  =  (1-a)T  ♦  aT,  where 


Tv  =  (P(pDfcv)  -wAv,  V|t_0 


Tv  =  (P(pD^v)  -WAv,  v 


1 1“0 


>  e  V, 

)  e  y. 
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Finally  da  not*  by  Y  the  sat  of  valuas  a  e  [0,1]  for  which  problem  (2.2)  la  solvable 


In  X  for  ovary  pair  <F,vo>  e  y.  Clearly  o  e  r,  because  for  this  value  of  the 
parameter  equation  (2.6)  because  the  linearised  Navier-Stokes  equation.  Let  us  verify 
that  r  la  open  and  closed. 

Y  la  ooen.  Let  a  e  Y  and  denote  by  G(F,v  )  =  v  the  solution  v  of  problea 

"  o  o 

T  v  -  (F,v  ).  Fro*  (2.3)  one  gets  G  e  L(y»X>  (2)  with  IGI  <  1.  Equation 

01  O  //A 

O 

T  .  v  -  (F,v  )  can  be  written  in  the  fora 
a  +c  o 

o  _ 

(2.7)  [1-e  G(T-T)]v  -  G(F,v  ). 

o 

Since  IG(T-T)I  v  <  IT-TI  ,  equation  (2.7)  is  solvable  for  |e|  <  It-TI*,1.,  (by  a 

X,  y  A»y  A,y 

Deusan  expansion). 

Y  la  closed.  Let  a  e  Y.  a  *  a  ,  and  let  V  be  the  solution  of 

m  mo  n 

Tv-  (F,v  ).  Fro*  (2.3)  one  has  <v  I  <  I (F,v  >1  .  Since  X  is  an  Hilbert  space 
a  n  o  n  a  o  y 

n  _ 

thsr*  exists  a  subsequence  ♦  v  e  X,  weakly  in  X#  Proa  T,  T  €  L(XiV)  one  has 

Tv  ♦  Tv,  Tv  ♦  Tv  weakly  in  V.  Hence  T  vy  ♦  T  v,  i.e.  Tq  v  -  (F,Vo). 
v  v  °v  o  o 

Let  us  now  return  to  problem  (1.1).  Define 

(2.8)  F(p,v)  =  P  {-p(v«V)v  +  X  (<v*V)Vp  *  (Vp*V)v)  + 

l2  1 

+  —  [(Vp*V)Vp-  f  (Vo*Vp)Vp  *  ApVp)  ♦  pf)  . 

P  P 

For  convenience  we  will  use  in  the  sequel  the  translation 

A 

(2.9)  p  -  p  +  O. 


Recall  that  P  is  a  given  constant.  To  solve  proble*  (1.1)  and  (1.2)  in  our 

2  2 

functional  fraaework  is  equivalent  to  finding  v  €  L  (0,TjD(A>),  v»  C  L  (0,T »•)  and 
C  e  L2(0,T»  H^>,  O'  e  L2(0,Ti  H1)  such  that 

%(PDtv)-|iAv  -  F(p,v), 

...  v.  -  v  (x), 

(2.10)  |t-0  o 

D^o  -X7o  -  -v*Vc, 

V-C  "  Vx)' 


(2)  Banach  space  of  linear  continuous  operator  fro*  V  into  X,  with  nor*  I  •  y  ^ 
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where  e  T  and  Oq(x)  -p  e  are  given.  Note  that  frost  the  above  conditions  on 


o  it  follows  that  a  e  c(0,Ti  «h). 

He  solve  (2.10)  by  considering  the  linearised  probla 
P(pD^v)-pAv  *  r(p,v)  =  F, 


(2.11) 


v|t-0  “  votx)' 

O^o  -  XAo  -  -v* Vo .  < 

°|t-0  ■  Vx)' 

and  by  proving  the  existence  of  a  fixed  point  (p,v)  ”  (p,v)  for  the  sap  (p,v)  ♦  (p,v) 
defined  by  (2.11). 

In  order  to  get  a  sufficiently  strong  estimate  for  the  linearised  equation  (211  «e 
take  in  account  the  particular  fore  of  the  date  v*Vcr.  As  for  estisate  (2.3)  we  wil. 
introduce  a  balance  parameter  e  >  0. 

Theoree  2.2  Assuee  that  v  e  L2(0,T»B^)  f>c(0,TiNg>  and  that 
a  e  L2(0,T»HjJ>  nc(O.Ti^).  Then  the  solution  o  e  L2(0,T»h*>,  o’  e  L2(0,T»H1)  of 
problee  (2. 11  )j.  (2. 11 )t  verifies  the  eatlsate 


(2.12) 


'°'c(o  t,  a2)  4  ,0,i2  3  ‘ 

c(o, ly  L2(0,T|H„) 


<  c,io  i2  +  c,  €~3t[ivi6  . 

2  O  2  w  —  |  j.  -  * 

C(0,tlH  ) 


+  IVol6 


1  1 
C( 0,T?H  ) 


♦  «^«2,  2  ]. 
I>  (0,T|H* ) 


♦  C  C  [Ivl2 

l  (0,T»H") 

for  every  positive  e  verifying 
(2.13) 

where  c,  is  the  constant  in  (2.16).  Here  ct,  c2,  c3  are  positive  constants  depending 
only  on  0. 

Proof.  The  existence  of  a  solution  0  in  the  required  space  follows  with  standard 
techniques  fro*  the  a  priori  bound  (2.12),  or  fro*  [9],  Vo X use  II,  Chapter  4,  Theoree  5.2, 


with  H  ■  H  .  Let  us  prove  ( 2. 12 ) s 
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By  application  of  the  operator  A  to  both  aides  of  (2.11)3,  by  multiplication  by 
Ac,  and  by  integration  over  Q  one  gets  V2  lAol2  +  (7(X&0-v*Vo) ,  VAo)  “  0.  Hence 

(2.14)  4  TT  140,2  *  XIVAol2  < 

2  at 

-  —  2— 

<  (iDvDol  ♦  IvO  ol )  IVAal  . 

16  " 

By  using  Sobolev'*  embeding  theorem  H  c>L  and  Holder's  inequality  it  follows  that 

(2.15)  4  S  IAcl 2+  XIVAol 2  < 

2  at 

_  1/  _  1/  _ 

<  CdDvr2  IDvl'2  |7ol  + 

+  lvl1  IVol’^2  1 74fl|1/2  )  I7A0I  . 

An  utilization  of  abc  <  (Be  3)a4  +  (e/2)  b2  +  (e/2)  c4,  e  >0,  leads  to 

(2.16)  4  -7  IAcl 2  +  XIVAol 2  <  c  elDvl 2  + 

2  at  1 

+  Cl  e!7Aol2  +  Ivl4  I  Vo  1 2  + 

+  c  e  IVAol2  *  Ivl2  IVol4  . 

Hence  for  e  verifying  (2.13)  one  has 

(2.17)  lAol2  ♦  X  I VAol 2  <  -f  (Ivl2  I  Vol 4  + 

at  e  3  1  1 

Ivl4  IVol2)  +  cge  (I Ovl 2  +  IVAol2), 

where  the  constants  c  depend  only  on  0.  'Biis  proves  inequality  (’. 12) .  Recall  that 
lolj  <  clAol,  lol3  <  cIVAol  . 

Remark .  One  could  also  consider  the  linearized  equation  0  +  v  Vo  -  XAo  -  0  instead 
of  (2.11)g>  then  estimate  (2.17)  holds  with  0  replaced  by  a  and  without  the  term 

—  2 

cge  I Viol  .  In  this  case  the  solution  0  of  the  linearized  problem  verifies  the  maximum 

principle  (which  doesn't  hold  for  the  solution  of  (2.11)j).  However,  the  linearization 

(2.11)3  seems  to  be  BK>re  in  keeping  with  the  linearization  (2.11  >^.  Besides,  the  maximum 

principle  will  be  recovered  for  the  solution  of  the  full  nonlinear  problem  (2.10)^. 

3.  The  nonlinear  problem.  local  existence.  He  will  not  take  care  of  the  explicit 

dependence  on  u,X,m,Mr  some  of  the  constants  c,Cy ,  depend  on  these  fixed  quantities.  In 

order  to  simplify  the  equations,  we  denote  by  XQ,  X^,  X^,...,  constants  depending  on  the 

norms  of  the  initial  data  lv  l_  and  lo  I,  . 

o  ▼  o  2 
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In  this  Met  ion  we  solve  (2.10)  by  proving  the  existence  of  a  fixed  point 
(P,v)  ■  (p,v)  for  system  (2.11).  Define 

. W, ' 

*  '*’">•"*>  •  ‘ !  c«  u.'i 1 

■s  ■  <; ,  olt.„  - 


«  2  c  In  I,,  io'l  ,  ,  <  K  , 

2  0  2  ^(O.T.H1)  ° 


i®  -  V  -  <  f  >  . 


Ctttj) 


where  C4  =  y[min  {u,M,(s»i2>/<4M2)}  ]  1  and  I  X/2c2  lo^  +  C  /«c2c4  IvqIv  IoqI2  • 
Here  c  -  c(fl)  is  a  positive  constant  such  that 

(3.1)  lv*wl,  <  c  lvlv  lw«2  ,  V  v  e  V,  w  e  H2  . 

Note  that  for  every  o  e  one  has  in 

(3.2)  2  <  p(t,x>  «  m  +  2  . 

Me  now  evaluate  the  I.2  norm  of  F  =  F(p,v).  By  using  Sobolev's  embeding  theorem 
16  " 

H  C» L  and  Holder's  inequality  one  easily  gets 


(3.3) 


lF(P,v)l2  <  C  Ivl2  I Dvl ^  C  Ivl2  I02ol  • 

•1D2oI1  +  C  Iv^  IDvl  1  I  DC  1 2  ♦ 

+  C  IDOl*  lD20l1  +  C  IDol®  + 

+  c  ifi2  . 


Consequently. 
(3.4)  lF(P,v)l2. 


L  (O.Tjl) 


<  C  (Iv  I  tlo  l,)7/2  TV2  ♦ 
O  V  o  « 


+  C  In  ll  T  ♦  c  ■  f  1 2 

°  L  (0,T>H) 


Hence,  by  using  (2.3),  it  follows  that  the  solution  v  of  (2.11).,,  (2.11)2  verifies 


(3<5>  'V'c(0,T,V>  +  ,V,22  2  *  ,V’,22  2 

C'  '  '  1  L(0,T| WT)  L  <0,T»O 


<  C,  lv  I2  ♦  K,(/T  +  T)  ♦  C  Ifl2 

4  o  V  1  . 2.-  _ 

L  ( 0 , T| 8)  • 

On  th«  other  hand,  from  (2.12), 

<3. 6)  lol2  ♦  lol2  < 

C(0,T>H*)  L  (O.TjHj,) 


<  {  C,  10  I2  ♦  K,  e'3T  ♦  K .  £  }  . 
4  O  2  4  4 


-1  2 

Now  we  fix  £  >  0  such  that  K  £  <  2  C,  lo  I, 

4  ^  0  2 


Finally,  by  fixing  a 

sufficiently  small  I  >  0,  it  follows  that  v  (  o  e  The  estimate  for  o 

follows  by  using  equation  (2.11 >3  and  (3.1).  The  estimate  for  the  sup  norm  of 


o  -  o  in  Q  is  proved  as  follows: 
o  t  . 

6  i/2 

Clearly,  lo(t)  -  o  I,  <  f  lo'(s)l,  ds  <  K  T  .  On  the  other  hand  lo(t)  -  o  I 
olio  o 

1/3  °  2/3  C((I) 

<  C.  lo(t)-o  I,  lo(t)-o  I,  ,  where  C5  depends  only  on  B»  recall  that 


5 
5/3 


o  1 
(0)  C*C(0). 


o  2 

Consequently 

1/3  1/6 


lo  -  a  I  <  C,  K 

°c(a>  5  ° 

T 


<C6  /2C2  ,0o'2  *  ,0o'2> 


2/3 


where  C,  «  C,((J)  is  a  positive  conatant,  such  that  lol  <  C.  lol,  ,  Vo  e  H  (fl). 

6  6  C(il)  6  2  " 

Hence,  by  choosing  (if  necessary)  a  smaller  value  for  T,  one  gets  lo  -  o  I  _  <  m/2  . 

Now  we  utilize  Shauder's  fixed  point  theorem.  Clearly  I  =  I,  ♦  lj  is  a  convex, 

2  2  2  2  —  — 
compact  set  in  L  (0,T>L  )  x  L  (0,T>L  ).  Let  us  denote  by  4  the  map  4(p,v)  •  <p,v). 


defined  by  (2.11).  Since  4  (K)  C  K,  it  is  sufficient  to  prove  that  4  i  K  ♦  K  is 

2  —  —  2  —  ~  2 
continuous  in  the  L  topology.  If  v  +  v  in  L  (Q),  p  ♦  P  in  L  (Q_),  is  follows  by 

n  r  n  * 

-  2  2  12 
compactness  arguments  that  vp  •*  v  weakly  in  L  (0,T|H  )  and  in  H  (0,T;L  ),  and  that 

-  -  2  2  12  — 

Vp  ♦  Vp  weakly  in  L  (0,TiR  )  and  in  H  (0,T;L  ).  In  particular  p  is  a  bounded 

n  V2te  2-e  3  n 

Z)  C»C0'a(QT)  '  ,  for  suitable  positive  Ej,  a. 


sequence  in  H 


(0,T:H 


( J )  a-Holder  continuous  functions  in  Q  . 
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-  -  -  -  -  y.  i 

Hence  p  ♦  p  uniformly  in  Q  .  Moreover,  v  and  Vp  are  bounded  in  H  *  (0,T|I  >  and 
v  t  n  n 

1/1  -  -  -  -  A 

in  HM0,T;H  )  respectively.  Hence  v  *  v  and  Vp  ♦  Vp  strongly  in  the  L 

n  n 

topology.  It  follows  from  (2.8)  that  r(p  ,  v  >  ♦  F(p,v)  as  a  distribution  in  Q-i 

n  n  * 

consequently  Ftp  ,  v  )  *  F<P,v)  weakly  in  L2(Q_),  because  P(p  ,v  )  is  a  bounded 
n  n  t  n  n 

sequence  in  this  space.  Analogously,  v^*  Vp^  ♦  v  •  Vp  strongly  in  L2 (Q^ > .  It  follows 

2  2 

frc*  the  linear  equations  (2.11)  that  vn  *  v  and  Pr  *  p  in  L  (Q^)  and  I.  (Q^ ) , 
respectively.  Hence  ♦  is  continuous.  This  finishes  the  proof  of  the  existence  of  a 
local  solution.  Uniqueness  will  be  proved  in  Section  S. 


a*  Global  solutions.  Asymptotic  behavior. 

In  this  section  the  constants  Cg  depend  at  most  on  Q  and  on  the  quantities 

A  A 

U,  1  and  p,  i.e.  on  the  total  amount  of  mass  |Q|  p.  We  assume  that 

(4.1)  lo  I,  <  (2C)"1  p, 

O  2 

hence  p/2  <  m  <  M  <  3  p/2.  Let  (P,v)  be  a  solution  of  (1.1).  From  (2.6)  for 
2  -1 

tQ  "  (4m  )  m  u,  and  from  (3.3)  one  gats 


(4.2) 


where  C 
one  gets 


—  ID  vl2  ♦  ^  Ivl2 
2  t  2  dt  V 

<  C  dvi2 

depends  only  on  fl ,  u ,  p  . 


+  ®-r  I AvI 2  < 

8  M 

♦  Id*  )  (lvl2  +  lalj)  +  C  Id j 
On  the  other  hand,  from  (2.17) 


+  C  Ifl2  , 
for  e  -  (2Cg)-1X, 


(4.3) 


lid2  +  ^  IVAd2  <  C  (Ivl®  +1  d®  ) 

at  2  1  2 


From  (4.2)  and  (4.3)  it  easily  follows  that 
~  (  If  'vl2  ♦  Idol2)  +  f  lDtvl2  + 


16M 


-  I  AvI 2  ♦  j  IVAol2  <  C(  I  vl®  +I6d®+lfl2) 
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(4.4) 


In  particular,  ainca  lAvl  >  Clvly  and  i  VAal  >  clAdl,  ona  haa 
OvlJ  ♦  lAol2)  <  -[C^-C^dvlj  ♦lAol2)2]. 

•  (Ivl2  +  lAol 2 )  +  C12  Ifl2. 

Hanca  (4.5 )1  below  holda  for  every  t  e  (0,+“[  if 


(4.5) 

In  fact,  if  C11  (  lv(t)l2  ♦  IAo(t)2!)2  -  cyo/2  *tt*t  b®'  ltro*  th,t 

“  (iv(t)l2  +  i Ao(t) I2)  <  0. 

Lat  ua  now  prova  the  laat  aasartion  in  thaoraa  A.  Under  the  hypothaaia  (4.5),,  it 
followa  froai  (4.4)  that 

Q 

(Ivl2  ♦  I Aol2 )  <  -  -5^  (Ivl2  ♦  IAOI2)  . 

thin  provaa  (1.6)  .  q 

5.  Onlguanaaa.  He  prova  that  the  aolution  (p,v)  of  problaa  (1.1)  ia  unique  in  the  claaa 
in  which  axiatanca  waa  proved;  aea  thaoraa  A.  Ha  raaark  that  aore  careful  calculationa 
lead  to  uniquenaae  in  a  larger  claaa. 

Lat  (P,v),  (p,v)  be  two  aolutione  of  problaa  (1.1),  (1.2)  and  put 
u  “  v-v,  n  -  p-p.  By  aubtracting  the  equationa  (2.10),  written  for  (p,v)  and  (p,v) 

raapactively,  and  by  taking  the  inner  product  with  u  in  H  one  gata 

J  (pu,u)  +  U  lul2  “  “  2  t»,vP«u2>  * 

♦  “  (Ap ,u2)  -  (w,Dtv*u)  ♦  (f-r,u)  . 

2  2 

By  uaing  (Ap,u  )  “  -(Vp.Du  >  ,  wa  a how  that 
(5.1)  “  —  (Pu,u)  ♦  ^  lul2  <  j  Ivl^lVpI^lul2  ♦ 

2 

X  2  2  —  2  2 

♦  ^  l?pl*  tul  ♦  ClDtvl  lul  ♦ 

j  i  Ant2  +  (r-r,u)  . 
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On  the  other  hand,  by  subtracting  aquations  (2.10)3  written  (or 
(P,v)  and  for  (p,v)  respectively,  and  by  taking  the  inner  product  with  An  in  L2 (Q) 
one  gets 

(5.2)  4  ,Vr'12  ♦  7  lAnl2  <  C  l?pl2  »ul2  ♦ 

2  at  2  m 

+  c  ivi2  I7ni2  . 

By  adding  (5.1)  and  (5.2)  it  follows  that 

(5.3)  l (pu,u)  ♦  IVtiI2J  ♦  y  lul2  ♦  |  lAnl2  < 

<  e^t)  (iui2  +  ivni2)  +  (p-r.u), 

where  9f(t)  ia  a  real  integrable  function  on  (0,T). 

On  the  other  hand,  by  using  Sobolsv's  sabeeding  theorems  and  Holder's  inequality 
2-12 

(and  also  ab  <  (a  M  b  )  the  reader  easily  verifies  that  given  e  >  0  there  exists  an 
integrable  real  function  82(t)  (dependant  on  P,  P,  v,  v  and  on  e)  such  that 

(5.4)  |(P-P,u)|  <  e2(t)  lul2  +  e  (lnl2+lul2  )  . 

2  -1 

By  using  lul  <  ■  (pu,u)  ,  (5.3)  and  (S.4)  it  follows  that 

^  [(pu,u)  +  I Vnl 2J  <  (Q^t)  +  fl2(t))  r(PU,u>  ♦  I Vri* 2 >  1  - 
Uniqueness  follows  now  fro*  Gronwall's  ltsu  and  from  U| t_Q  “  0,  “  0. 
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